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Abstract 

The solution of the Helmholtz equation in optical semiclassic ap- 
proximation is associated with the calculation of ray paths and ma- 
trices of variations. The transformation rules for elements of matrices 
on the boundaries of the waveguide are obtained. 

Key words: geometry acoustics, matrix of variations, boundary condi- 
tions 

Consider a two-dimensional underwater acoustic waveguide with the sound 
speed c = c (r, z) being a function of depth, z, and range, r, and curve bottom 
z = z b (r). 

It is well-known that in optical semiclassical approximation the sound 
wave field can be expressed through parameters of ray trajectories governed 
by the Hamilton equation [lj 
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with the Hamilton function 



H = — >y n 2 — p 2 . 



The variable p = n sini) is the ray pulse, •& is ray grazing angle, n = n (r, z) = 
Cq/c (r, z) is the refractive index, where c is a reference sound speed. 

To find a ray trajectory starting from the source point (ro, zq) at the 
angle $o we must solve the system (CQ) with the initial conditions z = zq, p = 



1 



n (r , Zq) sin^o- As there exist phase restriction for the variable z then we 
need in boundary conditions for variables in addition to the system ([I]). The 
conditions have the form 



p->p{l-2N g (N r t r /t g + N s )) 
z — >• z, 



where t = (t r , t z ) = (cos sin i?) is direction vector of the ray, N = (iV, 



D 1 * z j 



(cos a, sin a) is internal normal to the boundary at the point of ray reflection. 

Calculations of amplitudes of the sound pressure deals with the matrix 
of variations 

Qn Q12 



governed by the equation 
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Starting value for matrix q is the identity matrix. Coefficients of the lin- 
ear system ([2]) vary through the ray trajectory. When the ray reaches the 
waveguide boundary matrix q suffers a transformation of the form 



Q 



nq. 



The transformation matrix 



«11 
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where 



fi =f- 2iViV t = (t lr ,t l2 



(4) 



is the direction of a ray after reflection of it from the bottom, N t = (t, N) . 
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Figure 1: Coefficients k u (a), n 2 2 (b) and K\ 2 (c) for fixed ray gracing angles $ 
versus angle a; the curve number fc corresponds to the angle $ = — 90° + 10°/c; 
the curve at the point of reflection R = 0.02, derivatives n' z = 0.01, n' = 0. 



In particular case on the free surface z = and also on the horizontal 
bottom z = const we have 
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tz 
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As ft 12 7^ when ra' 2 ^ then horizontal planes scatter rays in inhomogeneous 
media. 

In homogeneous waveguide with n = const we obtain from ([3]) 

_s n +, + _2_ 
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To prove formula ([3]) let consider a narrow beam of rays, bounded by the 
rays with similar pulses p and p. Let the incident points of the rays be (r, z) 

and (r = r + 5r, z = z + 8z), direction vectors t and t, normals iV and N, 
respectively. 

From geometrical reason (Fig. l2b]) we receive 
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t r N z 
N, 



Az, 
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t r N r 
~ Az. 



First, we find the laws of reflection of the beam, neglecting the change 
of the normal vector at the site of incidence on the surface, that is, assum- 
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Figure 2: Scheme reflection of a narrow beam of rays from the curved surface 
with normal vector iV and the radius of curvature R. 



ing a flat surface. For the increment of z = z(r) — z(r) we have in this 
approximation 



Z\ (r + 5r) = Z\ (r) + — Zx5r — z (r) + — 5r, 

dr ti r 

d t z 
z\ (r + 8r) — z (r) + —z5r = z + -J-5r, 

Azi = Z\ (r + <5r) — 5 (r + <5r) = Az + ( — — -J- ] 5r 

\ti r t r J 

l + ^f^-^Az^-^Az. 

N t \ti r t r J J t\ r 
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For Ap = p (r) — p (r) we have analogously: 



Pi 



Pi (r + 5r) = pi (r) + —pidr = 

dr 

= p (r) - 2N Z (N r ^n 2 -p 2 (r) + p (r) N z ) + ^5r, 

V / ti r 

(r + 5r) = p (r + 5r) — 2N Z (N r ^/n 2 (r + 5r,z + Sz) — p 2 (r + 5r) + p (r + 5r) iV 2 



p (r + <5r) = p (r) + —pSr = p (r) H -5r, 

dr £ r 

Api = (r + 5r) — p (r + <5r) = 
~ Ap — 2N z N r ^a/ n 2 — p 2 (r) — a/w 2 (r + <5r, i + 5z) — p 2 (r + <5r) j 

-2iV 2 (p (r) — p Or + Sr)) N z + [ ^-p\ - -7-p ) 5r = 

\dr dr J 

Ap (l - 2N? + 2iV A |) + - | P (l - 2iVf + 2^iV r |)) fr 

+2iV 2 iV r — p= {n' r Sr + n' z Sz) = 
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Here we have used the identity 



1 - 2N 2 Z + 2N z N r — = - — , 

Now take into account the effect of the curvature of the surface at the 
point of reflection, neglecting variations in t and n. Let the radius of curva- 
ture of the bottom surface is equal to R, and she curvature & = 1/R. 

We have 



dN z = dsina = N r da, 



dN r = dcosa = —N z da. 
Differentiating (J4j), we obtain 
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-dt lz /2 = dN z (t,N\+ N z (t, w\ 



= dN z ^(t, + t z N z ^ + N z t r dN r = 
= (N r (2t z N z + t r N r ) - N 2 z t r ) da = -t lr da. 

Then (see Fig. [2b]) 



5r t 
da ~ -Ml = = -Sl-^rAz, 

N z N t 

so when multiplied by n we get 

Api = 2nt lr da ~ -2&n-^-Az. (7) 

After summation of (JHD, © and passing to the limit as Ap -> we 
obtain 
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